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By applying the recently developed nonperturbative functional renormalization group (FRG) 
approach, we study the interplay between ferromagnetism, quasi-long range order (QLRO) and 
criticality in the d-dimensional random field 0{N) model in the whole (A'', d) diagram. Even 
though the "dimensional reduction" property breaks down below some critical line, the topology of 
the phase diagram is found similar to that of the pure 0{N) model, with however no equivalent of the 
Kosterlitz-Thouless transition. In addition, we obtain that QLRO, namely a topologically ordered 
"Bragg glass" phase, is absent in the 3-dimensional random field XY model. The nonperturbative 
results are supplemented by a perturbative FRG analysis to two loops around d = 4. 



How the phase behavior and ordering transitions of a 
system are affected by the presence of a weak random 
field remains in part an unsettled problem. Heuristic 
and rigorous arguments show that the lower critical di- 
mension below which no long-range order is_possible is 
2 for the random field Ising model (RFIM) 1, 2] and 4 
for models with a continuous symmetry (RFO(iV)M with 

> 1) QIES' However, this leaves aside two questions: 
first, the nature of the critical behavior in random field 
models, a question connected to the breakdown of the 
so-called " dimensional reduction" (DR) property that re- 
lates the critical exponents of the RFO(A^)M to those of 
the pure 0{N) model in two dimensions less 0; and sec- 
ond, the possible occurence of a low-temperature phase 
with quasi-long range order (QLRO), i.e. , a phase char- 
acterized by no magnetization and a power-law decrease 
of the correlation fuctions, in models with a continuous 
symmetry Q . Progress has been made to better cir- 
cumscribe this latter point, ft has indeed been shown 
that QLRO is absent for N — 2 when disorder is strong 
and for > 3 for arbitrarily weak random field fs*!; but 
this still keeps open the cases iV = 2, 3 in the physical 
dimensions d = 2,3. 

Those questions are important because, on top of 
purely theoretical motivations, they concern the behavior 
of the known experimental realizations of random field 
models. This is the case for instance of vortex lattices 
in disordered type-H superconductors 0i 0|- such 
systems, the randomly pinned lattice of vortices can be 
mapped onto an " elastic glass model" 0, , whose sim- 
plest realization is the N ~ 2 RFXYM. The occurence 
of a phase with QLRO, termed "Bragg glass" has been 
predicted for the 3 — d version of the model |2| ■ Further 
theoretical support for this prediction has been given by 
a Monte Carlo simulation of the RFXYM 01 and by 
analyses of the energetics of dislocation loops 0, 0| . 

In this letter, we apply our recently developed nonper- 
turbative FRG approach of the RFO(A^)M fl3| to provide 
a unified picture of ferromagnetism, QLRO and critical- 
ity in the whole {N, d) diagram. We find that below a 
critical value Nr ~ 2.83.. and for d < A the model has 



a transition to a QLRO phase, both this phase and the 
transition being governed by zero-temperature nonana- 
lytic fixed points (FPs). The transition disappears below 
a lower critical dimension which we find around 3.8 for 
N = 2. Therefore, contrary to what is usually believed, 
no QLRO and no Bragg glass phase exist in the 3 — d 
RFXYM. We supplement our nonperturbative, but of 
course approximate, results by a perturbative FRG anal- 
ysis to two loops in d = 4 -|- e. The present approach 
allows us to discuss the DR property and its breakdown. 
We find in particular that the topology of the {N, d) 
phase diagram of the RFO(A^)M is similar to that of 
the pure 0{N) model, with however no equivalent of the 
{N = 2, d ^ 2) Kosterlitz-Thouless (KT) transition. 

Our starting point is the standard effective hamilto- 
nian for the RFO(iV)M in d dimensions with an A^- 
component field xi^) ^nd uncorrelated random fields 
taken from a Gaussian distribution with zero mean and 
variance A. After introducing n replicas in order to per- 
form the average over quenched disorder (taking at the 
end the limit n — s- 0), it can be rewritten as: 
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where we have introduced sources acting on each replica 
separately, which therefore explicitly break the permuta- 
tion symmetry between the replicas. 

We apply in this work the nonperturbative FRG for- 
malism for the RFO(A^)M recently proposed by us 0|. 
To keep the presentation short but sufficiently self- 
contained, we first sketch the main steps of the approach. 
It is based on an exact RG equation for the effective av- 
erage action rj;[{^a}]. interpolates between the bare 
action, eq.(E3), and the usual effective action (i.e., the 
Legendre transform of the partition function associated 
with eq.(^) as the running scale k moves from micro- 
scopic (fc — A) to macroscopic (fc 0) scale. It is built 
by integrating out fluctuations with momenta larger than 
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k and its flow obeys an exact equation, 

dtrk[{<f>a}] = iTr{ati?fe(g)(rf + IRkiq))-'} 

(2) 

where dt is a derivative with respect to t = ln(fc/A), Fj. ' 
is the tensor formed by the second functional derivatives 
of Tk with respect to the fields (pail), 1 is the unit ten- 
sor, and the trace involves an integration over momenta 
as well as a sum over replica indices and A^- vector compo- 
nents. Rk{q) is the infrared cutoff introduced to suppress 
the low-momentum modes. In practice, solving eg.© 
numerically requires the introduction of approximation 
schemes which amount to truncate the functional form 
of Ffe. Guided by the physics of the problem at hand 
and by the mounting work on the method [l^ . one can 
then formulate a nonperturbative RG description. For 
the RFO(iV)M we have argued that the two main ingre- 
dients allowing to study the long-distance physics are (i) 
the derivative expansion, which approximates the mo- 
mentum dependence of the (IPI) vertex functions 
and (ii) the expansion in increasing number of free replica 
sums, which is equivalent to include increasing-order cu- 
mulants of the renormalized distribution of the quenched 
disorder [l^ . 

The simplest nonperturbative FRG description of the 
RFO(A^)M relies on the following truncation: 



r, 



1 " 1 
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,.fc=i 



with one single wavefunction renormalization for all 
fields, Z^^is, evaluated at the field configurations that 
minimize the 1-repli ca p otential Uk (pseudo first-order 
derivative expansion the 1-replica Uk and 2-replica 

Vk parts of the effective potential are obtained from F^ 
taken for uniform fields. Inserting eq.© into eq.lO, leads 
to coupled partial differential equations for the func- 
tions Uk{p) and Vk{p, p' , z), where p = ^cjy^ and z = 
(p ■ 0'/-\/(4/5p'), and to a running anomalous dimension 
defined as ijk — —dt\ogZ,„^k- The sought FPs being at 
zero temperature, one also introduces a running tempera- 
ture Tk and the associated exponent 6k = dt logT^. This 
is most conveniently done by defining a renormalized dis- 
order strength A,„,fc = {2p^^k)~^d:,Vk{p,n,k, P^,k, z = 1), 
where p^^k corresponds to the minimum of Uk{p)- Defin- 
ing then Tk as Z'„^fcA:'^A/(AA„^fc), one can see that it 
reduces to the bare temperature T at the microscopic 
scale A. 

The flow equations can be expressed in a scaled form 
by introducing renormalized dimensionless quantities, 
Uk{p) = Tkk-'^Ukip), Vkip,p\z) = T^k-^Vk{prp',z), 
p = Z^^kTkk-'^'^-'^'ip: see eqs. (5,6) of ref. An ex- 

pression for rjk is derived by considering the flow of the 



transverse component of F^, 
and zero momentum: 



(2) 



evaluated for uniform fields 



rjk =^Cd/d{upp[vz - 4p„,Wpp']r7i3 ;^(2p„,Upp,0) 

+ 2u^pf.[m^^3(2pn,wpp, 0) + ml^{2p^upp, 0)]} 



(4) 



where C^^ = 2'^+^T:'^/'^T{d/2), Up 



and 



stand. 



respectively, for d^v and dpdpiv evaluated for fields 
that minimize Uk (the subscript k has been dropped for 
simplicity), and mt ^,(w,w') are dimensionless threshold 
functions described in [l3j. In eq.Q we have set for 
clarity Tk = T = 0, but the same FPs are reached when 
T > provided (as indeed found) that 6 = Ok^o > 0. As 
discussed previously , breakdown of DR is associated 
with the appearance of a strong enough nonanalycity in 
the renormalized 2-replica potential as the two replica 
fields become equal. This nonanalycity is related to the 
presence of metastable states in the renormalized random 
potential generated along the RG flow. 

We have numerically integrated the flow equations for 
Uk , Vk , and rjk for a variety of initial conditions and deter- 
mined in this way the FPs and the associated stability 
behavior (and critical exponents). To cover the whole 
(A^, d) diagram for continuous values of A^ and d with 
a tractable computational effort, we have used an addi- 
tional approximation that consists in expanding Uk and 
Vk around the field configuration p,„jij: Uk = U2(p— Pm.fc)^, 
Vk = voa{z) + vio{z){p + p' - 2p,„,fe) + V2a{z){p + p' - 
2pm,fe)^ + vo2{z){p — p'Y . One should point out that the 
present flow equations reproduce all 1-loop perturbative 
results in the appropriate region of the (A'^, c?) plane , in - 
cluding the FRG equation at first order in e = d — 4 . 
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FIG. 1: Predicted phase diagram for (a) the RFO(iV)M (non- 
perturbative FRG) and (b) pure 0{N) model (sketch). Re- 
gions I and II correspond to transitions to a ferromagnetic 
and a QLRO phase, respectively. 

The central results of our study are summarized on the 
phase diagram displayed in Fig.^. In region III there are 
no phase transitions, the RFO(A^)M always stays disor- 
dered. In region I there is a transition to a ferromagnetic 
phase at a critical point governed by a zero-temperature 
nonanalytic FP. The nonanalycity is strong enough (a 
"cusp") for breaking the DR predictions; however, when 
approaching a critical line Nc^Dnid) (not shown here, 
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but starting from N = 18 for d = 4 + e and going to 
d ~ 5 when TV = 10|), the critical exponents continu- 
ously tend to their DR value. Above the line, the stable 
(more precisely, once unstable) FP is nonanalytic, but the 
nonanalycity is now too weak to break DR. (FPs with a 
"cusp" can still be found, but they are unstable in sev- 
eral directions and correspond therefore to multicritical 
points, unreachable from generic initial conditions.) 

Finally, in region II we find two zero-temperature non- 
analytic FPs: one is attractive and describes a QLRO 
phase and the other one is once unstable and governs the 
transition between the disordered and the QLRO phases. 
The exponent rj characterizing the power-law decay of 
the connected correlation function is shown in Fig. El for 
a range of values of A^. One can see that for N less than 
a critical value Nc = 2.83... and for c? < 4, which cor- 
responds to region II, the two FPs coalesce for a value 
dic{N) which then determines the lower critical dimen- 
sion below which no phase transition is observed. The 
most striking outcome is that, contrary to what is usu- 
ally beheved, no QLRO, i.e., no Bragg glass phase, exists 
in d = 3 for the RFXYM: indeed, dic{N = 2) ~ 3.8. Our 
study of the RG flow for d = 3 and N = 2 shows signa- 
tures of a "ghost" nontrivial FP, presumably lying on 
the imaginary axis not too far from the physical plane 
of coupling constants (which could explain the behavior 
found in Monte Carlo simulation [loj. see Fig. EJ, but the 
flow goes at large distance to the trivial disordered FP. 
We shall come back to this point later. Note that an es- 
timate of the uncertainty of the present nonperturbative 
but approximate RG treatment is given by considering 
the point in Fig. ^ for d — 2: in this (probably most 
unfavorable) case, the theory predicts Nic ~ 1.15 instead 
of the exact result Nic = 1. As shown for the pure 0{N) 
modelfisj. this could be improved by solving the full first 
order of the derivative expansion. 




FIG. 2: Prediction for the d-dependence of the anomalous 
dimension rj for different values of (by steps of 0.2). Below 
Nc = 2.83... and for d < 4 two nontrivial FPs are found and 
coalesce for d = dic{N) (dots). 

It is instructive to contrast the behavior of the 
RFO(A)M with that of the pure model. The phase di- 
agram of this latter, derived from the solution of phe- 
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FIG. 3: Inverse correlation length vs disorder strength for the 
3-d RFXyM. The crossover around A = 5 is the remain of the 
FPs found above die. A parametrizes the initial conditions of 
the RG flow and 1/^ is estimated as in Ref. [T^ . 



nomenological RG equations (IS] and from known exact 
results, is shown in Fig. ^p. Regions I, II and II are the 
exact counterparts of those in Fig. QJl. Indeed, and al- 
though never mentioned, QLRO exists in the pure model 
even aside from the special KT transition for N — 2 
and d = 2. However, it occurs in an unphysical region, 
1<A<2, l<d<2! It remains that the topology of 
the RFO(Af)M phase diagram is very similar to that of 
the pure model. Again, there is no exact DR-like prop- 
erty, the lower critical dimension die (N) being shifted by 
an A^-dependent value between the two models. 

Another intriguing question raised by this similarity 
is the possible equivalence between the (Ac = 2.83..., 
d = 4) point of the RFO(Af)M and the {N^ = 2, d = 2) 
point of the pure model (compare Figs^ and b). Our 
nonperturbative FRG solution of the former shows a be- 
havior very reminiscent of the KT transition, but the 
approximation used is unable to rigorously locate a line 
of FPs as would occur for a KT transition. As studied in 
detail for the pure d = 2 XY model 0, ^3 ) ^ li^i^ of "al- 
most" FPs is found: this is enough to observe a behavior 
almost indistinguishable from the KT scenario, but not 
satisfying to prove the existence of a KT transition. To 
resolve this matter, we have considered the perturbative 
FRG analysis around d = 4. The 1-loop /3-function (or 
rather functional) for the renormalized 2-replica potential 
(to which our nonperturbative treatment exactly reduces 
near d = 4) can be shown to vanish identically in d = 4 
when A^ = Nc = 2.83... for arbitrary values of the dis- 
order strength, much like it does in the pure d = 2 XY 
model for arbitrary values of the temperature; but what 
about the 2-loop /3-function? To answer this question 
we have calculated (independently from the authors of 
the recent ref. ^) the /3-functions of the RFO(A^)M 
around d = 4 to 2 loops, starting with the nonlinear- 
sigma model associated with Eq. The calculation 

follows the treatment developed for the pure model [l9| . 
with the bare action expressed in terms of renormalized 
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dimensionless quantities, 



T •'^ a, 6=1 
h.d-2 r " y 



where the (A^— l)-component fields tTq (tTq = Hq) 
characterize the fluctuations in the directions orthog- 
onal to the external field h {h — zl^^ Z^^ H) while 
t = k'^^'^Z:j.^T and R{z) is the renor- 



terized by the propagator g''^^'^^ {q) 

Zl'"^ ZThlaaY^i^iLV- 



'^a = a/I - Z^-kI; 

malized 2-replica potential that we seek to relate to the 
bare term Ro{z) = Az via the counterterms S[R], such 
that Ro = k'^-'^{R + 6[R]). We then carry out the per- 
turbation expansion around the gaussian theory charac- 

tZrk^-'^iZ^q^ + 

-7T^Tr^)Sab where /i, = 1, . . .iV-1. 
The relevant interaction vertices are obtained by taking 
up to six functional derivatives of the action with respect 
to the TTa fields. They involve derivatives of R, i.e. , func- 
tions of ^. To handle the resulting functional diagrams 
we have followed the method developed for the FRG of 
disordered elastic systems Calculations are quite 

intricate and will be reproduced elsewhere. 

The final expression for the T = /?-function for 
R{z) can be put in the form dtR{z) = R'{l)(3[r] = 
R'{l){er{z) + R'{l)(3i[r] + R'{lfP2[r]), where R'{1) is the 
renormalized disorder strength and r(z) = R{z)/R'{l)- 
the 1-loop /?i [r] is given in refs [l^ l2l| and the 2-loop 
contribution reads 

2/32 H = yiyr'" - 3zr" - r'f{-yr" + zr' - I) - {N - 2) 
y2^//3 _ y(-3^^/ y _ 3)^//2 _ 2y(r' - z)r'r" + t/r'^ + 4r 



{N + 2)yr" - {3N - 2)zr' + 8K{N - 2)r 



where y = 1 — and a = lini(2^i-)(l — z'^^^'^r" (z). 
One can check that with the change of variable z — cos (f> 
and with K = 2"fa, the above equation is identical to the 
result recently obtained by Le Doussal and Wiese [l8j |. 
As expected, DR is recovered by setting the "anomalous" 
term a to zero. The coefficient K in the above expression 
is left so far as an unknown: it requires additional calcu- 
lations (in progress) in order to fully determine the 2-loop 
/3-functions. However, its value is not needed to answer 
the above question, i.e. , to know if for Nc — 2.83 . . . and 
d = 4 the /3-function for R{z) vanishes identically to two 
loops for arbitrary disorder strength R'{1). This latter 
property is only true if both f3i [r] and /32 [r] vanish iden- 
tically for the same function r{z). One can check that 
this is impossible, irrespective of the value of K. We can 
thus conclude that the special point {Nc — 2.83..., d — 4) 
does not correspond to a KT transition. Actually, as also 
obtained in our nonperturbative approach (see FigQJi), 



the only FP found for the (iVc, d = 4) model is at zero 
renormalized disorder strength, which allows a perturba- 
tive analysis of the two nonanalytic fixed points for d < 4 
and iV < TVc [13. 

One further consequence of the absence of a KT tran- 
sition is that, contrary to what occurs in the pure model 
near (TV = 2, d = 2), the line of lower critical dimension 
approaches the {N^ d — A) point with an infinite slope 
(compare Figs^ and b). This reinforces our finding that 
d = 3 is safely below the lower critical dimension of the 
RFXFM. Whether or not this conclusion implies that 
no Bragg glass phase should exist in disordered high- Tc 
superconductors is however an open question. The map- 
ping from these latter to the KFXYM is only valid in the 
low-T phase where an essentially elastic description can 
be used. It may well be that the mechanism by which 
QLRO is destroyed in the RFXFM, a mechanism that 
does not explicitly invokes the presence of dislocations, 
is specific to that model. 
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